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Abstract 

The neutrino mixing matrix is expanded in powers of a small parameter A, which 
\f~^ approximately equals to 0.1. The meaning of every order of the expansion is dis- 

cussed respectively, and the range of A is carefully calculated. We also present some 
qq ■ applications of this new parametrization, such as to the expression of the Jarlskog 

parameter J, in which the simplicities and advantages of this parametrization are 
^ \ shown. 
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5h ' In recent years, there have been abundant experimental data strongly suggest- 



ing the mixing of different generations of neutrinos, just analogous to that of 
quarks. The K2K [1] and Super-Kamiokande [2] experiments indicated that 
the atmospheric neutrino anomaly is due to the v ^ to v T oscillation with al- 
most the largest mixing angle of 6 atm m 45°. The KamLAND [3] and SNO [4] 
experiments told us that the solar neutrino deficit was caused by the oscilla- 
tion from v e to a mixture of and v T with a mixing angle approximately of 
9 sol ~ 34°. On the other hand, the non-observation of the V e to V e oscillation 
in the CHOOZ [5] experiment showed that the mixing angle 6 c h z is smaller 
than 3° at the best fit point [6,7]. 

These experiments not only confirmed the oscillations of neutrinos, but also 
measured the mass-squared differences of the neutrino mass eigenstates (the 
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allowed ranges at 3a) [6], 1.4xlO~ 3 eV 2 < Am 2 tm = |m|-m|| < 3.3xl0~ 3 e^ 2 , 
and 7.3 x 10~ 5 eV 2 < Am^, = (m 2 , — m 2 | < 9.1 x 10~ 5 ey 2 , where ± correspond 
to the normal and inverted schemes respectively. 

Like the Cabibbo-Kobayashi-Maskawa (CKM) [8,9] matrix for quark mixing, 
the neutrino mixing matrix is described by the unitary Maki-Nakawaga-Sakata 
(MNS) [10] matrix V, which links the neutrino flavor eigenstates u e , z/ M , v T to 
the mass eigenstates u±, z/ 2 , z/ 3 , 
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It is always feasible to parametrize the Majorana neutrino mixing matrix as 
a product of a Dirac neutrino mixing matrix (with three mixing angles and 
a CP-violating phase) and a diagonal phase matrix (with three phase angles, 
and only two of them are unremovable) [11]. In a form similar to the quark 
mixing matrix, the neutrino mixing matrix can also be written as follows 
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where Sj = sin^j, q = cos0j (for % = 1,2,3), 5 is the Dirac CP-violating 
phase and o"i, cr 2 , <r 3 are the Majorana CP- violating phases. If the neutrinos 
are of Dirac-type, the diagonal phase matrix on the right side hand of Eq. (2) 
can be rotated away by redefining the phases of the Dirac neutrino fields. 
The Dirac CP-violating phase is associated with the neutrino oscillations, CP 
and T violation. The Majorana CP-violating phases are associated with the 
neutrinoless double beta decay, and lepton-number- violating processes [12]. 

The three mixing angles 9 atm , 9 c h z , and 9 so [ are related to the three mixing 
angles #1, 82, and 63, which describe the mixing between 2nd and 3rd, 3rd 
and 1st, 1st and 2nd generations of neutrinos. To a good degree of accuracy, 
Oatm = Oii Ochz = O2, and 6 so i = O3. 

According to the results of the global analysis of the neutrino oscillation ex- 
periments, the elements of the modulus of the neutrino mixing matrix are 
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summarized as follows [13] 



0.70 



0.87 0.50 - 0.69 < 0.16 
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0.20 



0.61 0.34- 0.73 0.60 -0.80 



(3) 




0.63 0.36 - 0.74 0.58 -0.80 
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Quite different from the quark mixing matrix, almost all the non-diagonal 
elements of the neutrino mixing matrix are large, only with the exception of 
V e s. So it is unpractical to expand the matrix in powers of one of the non- 
diagonal elements, like the Wolfenstein parametrization of the quark mixing 
matrix [14]. Xing [15] has made the Wolfenstein-like parametrization for the 
neutrino mixing matrix, but they have to use much higher orders of the non- 
diagonal elements. In the quark mixing pattern, all the non-diagonal elements 
are small, so we may take it for granted that the mixing is a small modification 
to the unit matrix. But on the contrary, why could not we consider the large 
mixing as the common pattern, which is just the case in the neutrino mixing? 
So we may not expand the neutrino mixing matrix around the unit matrix. 

In this letter, we will just make an expansion of the neutrino mixing matrix 
based on the bi-maximal mixing pattern. Since there are two mixing angles 
near 45° {9\ « 45°, and #3 ~ 34°), the neutrino mixing matrix is not only the 
bi-large pattern as commonly said, but quite near the bi-maximal pattern, 
which reads 



Comparing with Eq. (3), we can make an expansion of V in powers of A, which 
satisfies 



where A measures the strength of the deviation of V e \ from the bi-maximal 
mixing pattern. Unlike the Wolfenstein parametrization of the quark mixing 
matrix, A here is at the diagonal element of the neutrino mixing matrix. Be- 
cause 0.70 < V e i < 0.87, A is a small positive parameter, which approximately 
equals to 0.1, and this expansion is reasonable and will converge quickly. Be- 
cause Q\ is quite near 45°, must be quite near v2/2 [1,2]. Then we can 
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set 

= V2/2 + aA 2 . (6) 



Also, since 82 is rather small (with the best fit point V e z = 0.045 [6]), we can 
set 

V e3 = b\\ (7) 



where a and b are both small parameters of order 1. 

Now we will calculate all the Sj and q (for i = 1, 2, 3) to the order of A 4 . From 
Eq. (7), s 2 = V e 3 = b\ 2 , we have 

1 



c 2 = y/1 - si = 1 - -tfX\ (8) 



From Eq. (6), we have 

sic 2 = = V2/2 + aA 2 , (9) 



using Eq. (8), we get 

Sl = ^1 + aA 2 + ^o 2 A 4 . (10) 
2 4 



Similarly 

c a = ^ - aA 2 - i(4v/2a 2 + v^A 4 , 
cs = f + A + ^ 2 A 4 , 

s 3 = — - A - V2X 2 - 2A 3 - -(12V2 + V2b 2 )X A . 
2 4 

(ii) 

Thus we obtain all the trigonometric functions of the three mixing angles. 

Now we can get all the elements of the neutrino mixing matrix straightfor- 
wardly, 
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fj,3 
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V2 



+ A, 



Vfc = ^ - A - V2\ 2 - 2A 3 - i(6v^ + V2b 2 )\\ 
K 3 = &A 2 , 



-1 + ^A + 1(2 + V2a - b)X 2 + 1(2^2 -2a- V2b)X 3 
+-(6 - 2V2a + 2a 2 - \/2a6 + b 2 )X\ 

= \ + ^fX - i(V2a + o)A 2 - I(2o - v^A 3 
-l(2a 2 - 2b + V2ab)\\ 



V2 



aA 5 



V Tl = ~ - - 1(2 - v^a + o)A 2 - l(2>/2 + 2a + v^A 3 
--(6 + 2V2a - V2ab)X 4 , 

V T2 = -\ - ^A - l(V2a + o)A 2 - \(2a - v / 2o)A 3 

+ 1(2& + V2ab-b 2 )X\ 
V T 3 = ^- aX 2 - \{2V2a 2 + V2b 2 )X 4 . 



(12) 
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Then we can expand the neutrino mixing matrix in powers of A, 
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(15) 



Now we will see the meaning of every order in the expansion of V. 

1. The term of A is the approximation of the lowest order, where the atmo- 
spheric and solar neutrino oscillations are both of the largest mixing angles of 



5 



45°. We call this the bi-maximal mixing pattern. 



2. The term of A 1 indicates the deviation of the neutrino mixing matrix from 
the bi-maximal mixing pattern. 

3. The term of A 2 shows the effect of the CP violation. Because the CP violation 
is described by the element V e3 [16], and in the terms of A and A 1 , V e3 = 0, 
the degree of the CP violation is of the order A 2 in our parametrization. 

4. The term of A 3 is the modification of higher order. 

So in our new parametrization the terms of A , A 1 , A 2 show the bi-maximal 
mixing pattern, the deviation from the bi-maximal mixing pattern, and the 
CP violation effect respectively. 

Next, we are going to determine the ranges of A, a and b. From the analysis 
above, we know that 



^ 2 + y^b 2 w 
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(16) 



The current experimental data of these three parameters are (the allowed 
ranges at 3a) [6,7] 

0.58 < sm6 atm <0.81, 
< sin 6^ < 0.16, 
0.48 < sin 9 sol < 0.61. (17) 

And the best fit points are [6,7] 



sin 8 atm = 0.72, 
sin chz = 0.045 
sin 6 so i = 0.55. 



From these three constraints, we can determine the ranges of the three pa- 
rameters A, a and b. First, we can determine the range of A. From Eq. (16), 
we have 
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Fig. 1. The allowed range of A, the solid curve is for 
0.65 - 2.83A - 4A 2 - 5.66A 3 - 12A 4 = 0, and the dashed curve is for 
0.28 - 2.83A - 4A 2 - 5.66A 3 - 12A 4 = 0. 

0.48 < — — A - V2X 2 - 2A 3 - -(12y/2 + y/2b 2 )X 4 < 0.61, 
2 4 

then 




0.28 - 2.83A-4A 2 -5.66A 3 -12A 4 , 2 0.65 - 2.83A - 4A 2 - 5.66A 3 - 12A 4 



<b z < 



A 4 A 4 
From the right part of the inequality, we have 



0.65 - 2.83A - 4A 2 - 5.66A 3 - 12A 4 > 0. (19) 

The allowed range of A is shown in Fig. 1, and we can see that —0.67 < A < 
0.17. 

From the left part of the inequality, if 



0.28 - 2.83A - 4A 2 - 5.66A 3 - 12A 4 > 0, 

we have —0.62 < A < 0.08, which does not agree with the value of V e \. So we 
must set 



0.28 - 2.83A - 4A 2 - 5.66A 3 - 12A 4 < 0, (20) 

and thus the inequality satisfies automatically. The allowed range of A is shown 
in Fig. 1, and we can see that A < —0.62 or A > 0.08. Summarizing these re- 
sults, we get 0.08 < A < 0.17, which is consistent with the primary estimation 
in Eq. (5). So it is reasonable and practical to expand the neutrino mixing 
matrix in powers of A. 
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Fig. 2. The range of b. 
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Fig. 3. The range of a. 



Now we can determine the range of b. Because bX 2 = sin# 2 = sin 9 c h z , using 
Eqs. (17), (18) and 0.08 < A < 0.17, we have that b = 0.045/A 2 . The range of 
b is shown in Fig. 2, and we can see that 1.56 < b < 7.03. 



Eq. (18), we have 0.58 < Sic 2 = siyl — s 2 < 0.81, with the best fit point 

0. 72. Thus y/2/2 + aX 2 = 0.72, so a = 0.01/A 2 . The range of a is shown in 
Fig. 3, and we can see that 0.35 < a < 1.56. 

In our new parametrization, several other corresponding observable quantities 
associated with the elements of the neutrino mixing matrix can be expressed 
in relatively simple forms. From the ranges of A, a and b, we can determine 
the ranges of these observable quantities. 

1. The Jarlskog parameter J [17]. J is the rephasing-invariant measurement 
of the lepton CP violation. The Majorana CP-violating phases can be re- 
moved away by redefining the phases of the Dirac fields, so only 5 is associ- 
ated with the CP violation. J = Im(V e 2V fl3 V* 3 V* 2 ) = sis 2 s 3 cic 2 c 3 sin 5. In our 
parametrization, J can be expressed in a very simple form (to the order of A 4 ) 



Similarly, in the case of a, 



>/2/2 + aX 2 



S\C2- Using Eq. (17) and 





(21) 
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Fig. 4. The range of J. 

Because si, S3, c\, C3 all have the factor y/2/2, there are four \/2/2 in J. 
So the degree of the lepton CP violation is suppressed four times. Again, J 
is suppressed by the factor b\ 2 = 0.045 [6]. Using 0.08 < A < 0.17, 1.56 < 
b < 7.03, we can determine the range of J in Fig. 4, and we can see that 
J « 0.00996 ~ 0.01096 (here we take sin 6 ~ 1). 

2. The effective Majorana mass term (m) ee . In the neutrinoless double beta 
decay, the effective Majorana mass term is defined as follows 

(m) ee = K^e 2 " 1 + m 2 V 2 2 e 2i ^ +m 3 V£e Mos |, 

where <7i, cx 2 , 03 are the Majorana CP-violating phases [12]. Using Eq. (12), 
we get 

(m) ee = |^(m ie 2iCTl + m 2 e 2iCT2 ) + (v^A + A 2 )(m ie 2iffl - m 2 e 2iCT2 ) 

+6 2 A 4 (m 3 e 2iCT3 - m 2 e 2ia2 )\. (22) 

We can see that the coefficients of the three terms show the influences of the 
three orders of A. Only m\ and m 2 are important to the value of (m) ee , and 
influence of 7723 almost vanish if the masses of the three mass eigenstates are 
nearly degenerated, because the coefficient 6 2 A 4 is of 10 -3 . 

3. The effective mass terms of neutrinos. The effective mass terms of neutrinos 
can be defined as follows (here we take electron neutrino for example.) 

(m) 2 ee m\\V e i\ 2 + m 2 2 \V e2 \ 2 + m 2 3 \V e3 \ 2 . 
Using Eq. (12), we get 

( m )l = h m \ + mj) - (y/2\ + A 2 )(m 2 - m\) + 6 2 A 4 (m 2 - m 2 ). (23) 
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Again, the coefficients of the three terms show the influences of the three 
orders of A. Noting that Am^ = \m\ — m\\ and Am^ ffl = |m| — we can 
rewrite Eq. (21) into 

(m)l = ml ± [i - (V2A + A 2 )]Am^, ± & 2 A 4 Ar^ m . (24) 

We can see from Eq. (24) that (m)\ is directly related with the masses and 
the mass-squared differences of neutrinos. So these two kinds of different ob- 
servable quantities are associated together in our parametrization. If we can 
separately measure Am^ m , Am^, and (m) 2 e to a good degree of accuracy, 
we can fix the value of mi, which will help us determine the absolute mass of 
neutrino ultimately. 

In summary, although all kinds of parametrization of the neutrino mixing ma- 
trix are mathematically equivalent, and applying any of them does not have 
any specific physical significance, however, it is quite likely that some par- 
ticular parametrization does have its usefulness and advantages in analysis 
of various experimental data. Furthermore, we can express other observable 
quantities in a simple and transparent way, and can link several different kinds 
of observable quantities together. This is the purpose of our new parametriza- 
tion, and we hope that this new parametrization will be useful in the phe- 
nomenology of neutrino physics. 
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